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Using a sample of 68.3 million Kl -k^-k^-k^ decays collected in 1996-1999 by the KTeV (E832) 
experiment at Fermilab, we present a detailed study of the Kl Tr^Tr'V Dalitz plot density. We 
report the first observation of interference from Kl tt'^tv'tt^'' decays in which tt'^-k~ rescatters 
to 7r°7r° in a final-state interaction. This rescattering effect is described by the Cabibbo-Isidori 
model, and it depends on the difference in pion scattering lengths between the isospin 7 = and 
7 = 2 states, ag — 02- Using the Cabibbo-Isidori model, and fixing [ag — a2)m^+ = 0.268 ± 0.017 as 
measured by the CERN-NA48 collaboration, we present the first measurement of the Kl — > 7r''7r''7r'' 



quadratic slope parameter that accounts for the rescattering effect: 
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(-hO.59 ± 0.20»tat ± 

0.48 syst ± 1.06ea;t) X 10 ^, where the uncertainties are from data statistics, KTeV systematic 
errors, and external systematic errors. Fitting for both ^ooo and ao — 02, we find /looo = (—2.09 ± 
0.62,t,t ± Q.72syst ± 0.28e:rO X 10-^ and (ao - 02) = 0.215 ± 0.014,tat ± 0.025,^.* ± 0.006,^ ; 
our value for ao — a2 is consistent with that from NA48. 

PACS numbers: 13.25.Es, 14.40.Aq 
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I. INTRODUCTION 

The amphtude for the Kl — > Tr^Tr^Tr" decay includes 
contributions from two sources. The first source is from 
intrinsic dynamics that represent a decaying directly 
into the Tr^Tr^Tr" final state. The second contribution is 
from the decay Kj^ — 7r+7r~7r° followed by a rescattering, 
7r"'"7r~ tt'^tt'^. The amplitudes from these two contri- 
butions result in a small (~ 1%) interference pattern in 
the Dalitz plot density. 

The Dalitz plot density corresponding to the intrin- 
sic K —f Zt: decay amplitude is approximately described 

by Q 

(1) 

where 



Xd = (5i-S2)/(\/3m2 + ) (2) 

Yd = (s3-so)/™^ (3) 

= {Pk-P^)\ i = l,2,3 (4) 

So = (si + S2 + S3)/3 , (5) 



and Pk and Pi are the four-momenta of the parent kaon 
and the three pions. The linear g,j parameters and the 
quadratic /i, k parameters are determined experimentally. 
For the specific case of Kl tt^tt^tt'^ decays, the linear 
terms vanish, h = 3k and the intrinsic Dalitz plot density 
reduces to 0, Q 

\M^^o\^ ^ 1 + hoooRl , (6) 

where 

Rl = Xl+ Yl (7) 

= [sq - (S1S2 -f S1S3 -I- S2S3)/3] , (8) 

and /looo is the quadratic slope parameter. With |/iooo| < 
10^^ and P^Y) < 2.5 from kinematic constraints, the vari- 
ation of lA^oooP over the entire Dalitz plot is less than 
2%. 

There are two previous measurements of /iqoo- The 
first reported measurement, from Fermilab experiment 
E731 U, is based on 5 million recorded Kj^ t^'^tPtP 
decays with an i?||-resolution (j{P?d) ^ 0.1 deter- 
mined from Monte Carlo simulations. CERN experiment 
NA48 [| used nearly 15 million decays with a{Rjj) ~ 
0.03. The average of these two results is /iooo(PDG06) = 
(—5.0 ± 1.4) X 10"'^ [l|. Here we report a more precise 
result from KTeV based on 68.3 million decays with 
a{Rjj) ~ 0.014. Compared with the previous measure- 
ments of ft-ooo, a major difference in the KTeV anal- 
ysis is that we take into account the contribution from 
Kl — *■ 7r"'"7r~7r'^ decays in which 7r~'"7r~ — » tt'^tt'^ in a final- 
state interaction. 

A full treatment of rescattering in X — > 37r decays, in- 
cluding higher order loop corrections, is given by Cabibbo 



and Isidori This model, referred to hereafter as 

CI3PI, describes a "cusp" in the region where the min- 
imum tt'^tt'^ mass is very near 2m^+ : a cusp refers to 
a localized region of the Dalitz plot where the density 
changes very rapidly. A precisely measured shape of this 
cusp can be used to measure the difference in pion scat- 
tering lengths between the isospin 7 = and 1 = 2 states, 
ao-a2. In 2006, the CERN-NA48 collaboration reported 
the first observation of a cusp in — > tt^tt^tt^ decays. 
The interference effect is from the decay n^n^Tr^ 
followed by rescattering: tt+tt" — > tt^tt". They reported 
(ao — a2)m^+ = 0.268 ± 0.017 Q, in excellent agreement 
with the prediction of Chiral Perturbation Theory [l, [§] . 

Compared to ■n^ir^iT^ decays, a much smaller 

cusp is expected in Kl Tr^Tr^Tr*^ decays, and here we 
report the first such observation as part of our measure- 
ment of the quadratic slope parameter. The expected 
distortion of the Kl tt^tt^tt" Dalitz plot is shown in 
Fig. [1^ using /iqoo = —0.005 and no contribution from 
rescattering, and in Fig. [TJj using /looo = and CI3PI @ 
to model rescattering from 7r+7r~7r°. 

The effects of rescattering and a negative value of ft-ooo 
both result in the Dalitz plot density dropping slowly 
as increases. The maximum variation is only a few 
percent. The main feature that separates these two 
effects is that while the quadratic slope parameter re- 
sults in a smooth linear function of , the rescattering 
from Kl 7r+7r^7r° results in a much sharper fall-off 
near some of the Dalitz plot edges (see "cusp" labels in 
Fig- [It)- Also note that the Dalitz plot density in Fig.[T^ 
is azimuthally symmetric, while the density in Fig. [T|3 is 
azimuthally asymmetric. This cusp will become more 
apparent when we examine the minimum tt^tt^ mass in 

§IYI1 

The outline of this report is as follows. The KTeV 
detector and simulation are described in § [TIllIIII The 
reconstruction of Kl — > tt^tt^tt^ decays and the deter- 
mination of the Dalitz plot variables {Xd,Yo) are pre- 
sented in § [IVl § |V] describes the fitting technique used 
to extract the quadratic slope parameter (/iqoo) and the 
difference in scattering lengths (oq — 02). Systematic un- 
certainties are described in § IVIIl and § IVIIII presents 
results for /looo with ao — 02 fixed to the value measured 
by the NA48 collaboration. In § IIXI both the quadratic 
slope parameter and the difference in scattering lengths 
are determined simultaneously in a two-parameter fit of 
the Kl 7r°7r°7r° phase space. 



II. EXPERIMENTAL APPARATUS 

The KTeV detector has been described in detail 
elsewhere [l3, [HI- Here we give a brief description of 
the essential detector components. An 800 GeV proton 
beam incident on a beryllium-oxide (primary) target pro- 
duces neutral kaons along with other charged and neu- 
tral particles. Sweeping magnets remove charged parti- 
cles from the beamline. Beryllium absorbers 20 meters 
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FIG. 1: Expected deviation from Kl Tr'^Tr^Tr" phase-space based on (a) ftooo = —0.005 and no contribution from rescattering, 
and (b) Kl — > 7r^7r~7r° with rescattering as calculated by Cabibbo and Isidori ^6|, with feooo ~ 0. The intensity scales are 
slightly different to better illustrate the relative shapes. 



downstream of the target attenuate the beam in a manner 
that increases the kaon-to-neutron ratio. A coUimation 
system results in two parallel neutral beams beginning 
90 meters from the primary target; each beam consists of 
roughly equal numbers of kaons and neutrons. The fidu- 
cial decay region is 121-158 meters from the target, and 
the vacuum region extends from 20-159 meters. A regen- 
erator, designed to produce Ks decays for the Re{t' /e) 
measurement, alternates between the two beams. The 
other neutral beam is called the vacuum {Kl) beam. 
Only Kl tt'^tt'^tt'^ decays from the vacuum beam are 
used in this analysis. 

The KTeV detector (Fig. ^ is located downstream 
of the decay region. The main element used in the anal- 
ysis is an electromagnetic calorimeter made of 3100 pure 
cesium iodide (Csl) crystals (Fig. For photons and 
electrons, the energy resolution is better than 1% and the 
position resolution is about 1 mm. The Csl calorimeter 
has two holes to allow the neutral beams to pass through 
without interacting. 

A spectrometer consisting of four drift chambers, two 
upstream and two downstream of a dipole magnet, mea- 
sures the momentum of charged particles; the resolution 
is C7p/p ~ [1.7 ® (p/14)] X 10~^, where p is the track mo- 
mentum in GeV/c. Bags filled with helium are placed 
between the drift chambers and inside the magnet, re- 
placing about 25 meters of air. A scintillator "trigger" 
hodoscope just upstream of the Csl is used to trigger 
on decays with charged particles in the final state. The 
KTeV beamline has very little material upstream of the 
Csl calorimeter, thereby reducing the impact of external 
photon conversions (7X Xe'^e~). The total amount 



of material is 0.043 radiation lengths, about half of which 
is in the trigger hodoscope. Eight photon-veto detectors 
along the decay region and spectrometer reject events 
with escaping particles. 

An electronic trigger for Kl — > tt^tt^tt" decays re- 
quires at least 25 GeV total energy deposit in the Csl 
calorimeter, as well as six isolated clusters with energy 
above 1 GeV. For Kl tt^tt^tt^ decays that satisfy the 
energy and vertex requirements (§ IIV|) . approximately 
9% of these decays satisfy the six-cluster trigger, and 
20% of the six-cluster events were recorded for analy- 
sis. The combined six-cluster data from three run peri- 
ods (1996, 1997, 1999) has nearly 400 million recorded 
events. Kl 7r°7r°7r° decays are recorded for use as a 
high-statistics crosscheck on the Monte Carlo (see below) 
determination of the acceptance in the i?e(e'/e) analysis 
(10] . This sample is ideal to study the Dalitz density. 



III. MONTE CARLO SIMULATION 

A Monte Carlo simulation (MC) is used to determine 
the expected Kl — > tt^tt^tt^ Dalitz plot density that 
would be observed without the contribution from rescat- 
tering and with /igoo = 0; i.e, pure phase space. The 
Kl TT^TT^TT^ dynamics are determined from devia- 
tions between the observed Dalitz density (Fig. O and 
that from the phase-space MC. The simulated phase- 
space density accounts for detector geometry, detector 
response, trigger, and selection requirements in the anal- 
ysis. 



The simulation of Kl tt tt tt decays begins by se- 
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lecting the kaon momentum from a distribution measured 
with K — > 7r"'"7r~ decays. Each simulated kaon undergoes 
scattering in the berylhum absorbers near the target, and 
kaons that hit the edge of any coUimator are either scat- 
tered or absorbed. For kaons that scatter from a col- 
hmator edge, the Kl-Ks mixture has been determined 
from a study of — > 7r~'"7r~7r'^ and tt+tt" decays. 

After generating a kaon trajectory downstream of the col- 
hmator, each photon from Kl — > tt^tt^tt^ — > 67 is traced 
through the detector, allowing for external (7 e+e~) 
conversions. The secondary electron-positron pairs are 
traced through the detector and include the effects of 
multiple scattering, energy loss from ionization, and 
bremsstrahlung. The effects of accidental activity are in- 
cluded by overlaying events from a trigger that recorded 
random activity in the detector that is proportional to 
the instantaneous intensity of the proton beam. 

For photons and electrons that hit the Csl calorimeter, 
the energy response is taken from a shower library gen- 
erated with GEANT each library entry contains the 
energy deposits in a 13 x 13 grid of crystals centered on 
the crystal struck by the incident particle. The shower 
library is binned in incident energy, position within a 
crystal, and angle. 

For both data and MC, the energy calibration for 
the Csl is performed with momentum-analyzed electrons 
from Kl — > TT^e^i' decays. To match Csl energy resolu- 
tions for data and MC, an additional 0.3% fluctuation is 
added to the MC energy response. The data also show 
a low-side response tail that is not present in the MC, 
and is probably due to photo-nuclear interactions in the 
Csl calorimeter. As explained in Appendix B of 
this energy-loss tail has been accurately measured with 
electrons from Kl — * TT^e^i/ decays, and this tail is em- 
pirically modeled in the simulation with the assumption 
that the energy-loss tail is the same for photons and elec- 
trons. Losses up to 40% of the incident photon/electron 
energy are included in the model. 

The Csl position resolution is measured with precise 
electron trajectories in Kl n^e^v decays. The po- 
sition resolution for the MC is found to be nearly 10% 
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FIG. 2: Plan view of the KTeV (E832) detector. The 
evacuated decay volume ends with a thin vacuum window at 
Z = 159 m. Only decays from the vacuum beam are used to 
measure /looo. 
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FIG. 3: Layout of Csl calorimeter. The two neutral beams 
go through the 15 x 15 cm^ beam holes (into page) shown by 
the two inner squares. The fiducial cut, indicated by the dark 
lines, rejects Kl — *■ tt^tt^tt'' decays in which any photon hits 
a crystal adjacent to a beam hole or at the outer boundary. 



worse than for data, requiring that the MC cluster po- 
sitions be "un-smeared" to match the data resolution. 
The un-smearing is done for each simulated photon clus- 
ter by moving the reconstructed position closer to the 
true (generated) position in the Csl calorimeter. The 
position un-smearing fraction is 0.07 -I- 0.0016i?-y, where 
is the photon energy in GeV. 
Nearly five billion Kl tt'^tt^tt^ decays were gener- 
ated by our Monte Carlo simulation. More than 90% 
of the generated decays are rejected by the geometric 
requirement that all six photons hit the Csl calorime- 
ter; about 2/3 of these six-cluster events are rejected by 
the selection requirements described below in § IIVI The 
resulting sample of 124.9 million reconstructed Kl 
7r°7r°7r" decays corresponds to 1.8 x the data statistics. 



IV. RECONSTRUCTION OF Kl TrV^Tr" 
DECAYS 

The reconstruction of Kl — > vr^Tr'^Tr'^ — > 67 is based on 
measured energies and positions of photons that hit the 
Csl calorimeter. Exactly six clusters, each with a trans- 
verse profile consistent with a photon, are required. The 
cluster positions must be separated by at least 7.5 cm, 
and each cluster energy must be greater than 3 GeV. 
For the two nearest photon clusters in the Csl calorime- 
ter, we require that the minimum-to-maximum photon 
energy ratio is greater than 20%; this requirement elim- 
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inates the most extreme cases of overlapping clusters in 
which an energetic photon lands very close to a photon 
of much lower energy. The fiducial volume is defined by 
cluster positions measured in the calorimeter; we reject 
events in which any reconstructed photon position is in a 
crystal adjacent to a beam-hole or in the outermost layer 
of crystals (Fig. [3]). 

To remove events in which the kaon has scattered in the 
collimator or regenerator, we define the center-of-energy 
of the six photon clusters to be 

Xce = ^ XiEj ^ Ei Uce = ^ ViEi/ ^ Ei (9) 

i i i i 

where Xi, yi are the measured photon positions in the Csl 
calorimeter, Ei are the measured photon energies, and 
the index i — 1,6. The coordinate XccUce is the point 
where the kaon would have intercepted the plane of the 
calorimeter if it had not decayed. The size of each beam 
at the Csl calorimeter is about 10 x 10 cm^; the center- 
of-energy, measured with ~ 1 millimeter resolution, is 
required to lie within an 11 x 11 cm^ square centered on 
the kaon beam. 

Photons are paired to reconstruct three neutral pions 
consistent with a single decay vertex. There are 15 pos- 
sible photon pairings for a K]^ — Tr^Tr^Tr" decay. To 
select the best tt^tt^tt^ pairing, we introduce a "pairing- 
X^" variable {xlair)^ which quantifies the consistency of 
the three n° vertices. To ensure a reliable reconstruction 
of the Dalitz variables, we require that the smallest of 
the 15 Xpair values is less than 10 (the mean Xpair 3), 
and also that the second smallest Xpair value is greater 
than 30. The location of the kaon decay vertex (z) is 
determined from a weighted average of the tt*^ vertices. 

The main kinematic requirement is that the invari- 
ant mass of the 7r°7r°7r° final state be between 0.494 
and 0.501 GeV/c^, or nearly zLAa. Figure [4^ shows the 
7r°7r°7r"-mass distribution for data and the MC. The mass 
side-bands are due to Kl — > tt^tt^tt^ decays in which 
the wrong photon pairing is found in the reconstruction. 
The fraction of reconstructed events outside the invari- 
ant mass cut is 0.21% for data and 0.20% for the MC, 
confirming that the MC provides an excellent description 
of the data. 

Additional selection requirements are that the energy- 
sum of the six photon clusters lie between 40 and 
160 GeV, and that the reconstructed decay vertex is 
within 121-158 meters from the primary target. To pre- 
vent an accidental cluster from faking a photon, we use 
the energy-vs-time profiles recorded by the Csl calorime- 
ter. For each photon candidate, the Csl cluster energy 
deposited in a 19 nanosecond window before the event 
must be consistent with pedestal. To limit the effect 
of external photon conversions in the detector material 
(7X — > Xe^e^), we allow no more than one hit in the 
scintillator hodoscope that lies 2 meters upstream of the 
Csl calorimeter. 

To improve the resolution of the Dalitz plot parameters 
{X £) , Yd , Rjj) , the cluster energies are adjusted for each 
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FIG. 4: (a) Invariant Tr^Tr^Tr" mass with all selection require- 
ments except for the 7r'^7r''7r''-mass and Xe- The Tr^Tr'^Tr" mass 
resolution (from Gaussian fit) is 0.94 MeV/c^. (b) shows Xe 
distribution with all other selection requirements; last bin in- 
cludes all events with Xe > 200. Dots are data and the 
histogram is MC. Vertical arrows show the selection require- 
ments. 

event by imposing kinematic constraints to minimize 

i=l '<■ 

where Ei are the reconstructed cluster energies, ct^ are 
the energy resolutions, and i?/** are the six fitted clus- 
ter energies. The impact of cluster position resolution 
on the Dalitz parameters is much smaller than that of 
the energy resolution, and therefore the cluster positions 
are fixed in the minimization. The four kinematic con- 
straints are M(,~f = Mk and M^^ = TOttO for each of the 
three neutral pions. With four constraints and six un- 
knowns in Eq. I10[ the minimization has two degrees of 
freedom. Events with x|. < 50 are selected for the analy- 
sis. The minimization of x|; improves the i?|)-resolution 
from 0.070 to 0.014. Fig. shows a data-MC compar- 
ison of the xl; distribution. The fraction of events re- 
moved by the Xe cut is 0.43% for data and 0.47% for 
the MC; this slight disagreement will be addressed in the 
evaluation of systematic uncertainties. 

After all reconstruction and selection requirements 
there are 68.3 million tt^tt^tt" events. The two- 

dimensional Dalitz plot distribution for this sample is 
shown in Fig. [5] with no acceptance correction, and pro- 
jections onto and the minimum tt^tt^ mass are shown 
in Fig. [6K,b. The ~ 5% variation across the Dalitz plot 
is mainly from the detector acceptance; this variation is 
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is described within the framework of a single-parameter 
fit for hoQQ with the value of oq — a2 fixed by an external 
measurement (from NA48). However, this fitting tech- 
nique works the same way when both /looo and ag — 02 
are floated in the fit. 



TABLE I: Parameters and their values used in the CI3PI 



model. Kstt refers to Kl — * 
to isospin I — and 1 = 2. 



n n TT , subscripts and 2 refer 



parameter 



value 



Kstt linear slope {g+-o) 
Kstt quadratic slope [h^ — 0) 
aom^+ at it'^h~ thresh 
(ao — a2)m^-ir at tt+tt" thresh 
Effective ranges (ro,r2) 

Isospin breaking parameter (e) 



0.678 ± 0.008 fl] 
0.076 ± 0.006 [IJ 
0.216 ±0.013 [13] 
0.268 ±0.017 [7] 
1.25 ±0.04, 1.81 ±0.05 [6] 
0.28 ±0.03 [6] 
0.065 [7] 



FIG. 5: Dalitz plot density. Yd vs. Xd, for 68.3 million 
Kl — * TT^vr'^Tr'' decays in the KTeV data sample after all 
selection requirements. The color-scale at right shows the 
number of events in each 0.05 x 0.05 pixel. The reconstruction 
resolution on Xd and Yd is cr '--^ 0.01 as determined by the 
MC; the box in the lower-left corner shows 10a x lOcr for 
illustration. 



nearly an order of magnitude larger than the variations 
from rescattering (C13PI) and the quadratic slope pa- 
rameter. Also note that the uncorrected phase space dis- 
tribution has a minimum at the center of the Dalitz plot, 
while the expected effects from physics (Fig. [T]) result in 
a maximum at the center. An accurate —> tt^tt'^tt'^ 
simulation is therefore critical to this measurement. 



V. FIT FOR hooo AND ao 



02 



In the previous measurements of the quadratic slope 
parameter /iqoo [1, the i?|)-distribution for data 
(Fig. [5^) was compared to the R 



I) distribution from a 
simulated Kl tt^tt^tt^ sample with hooo = 0. Nor- 



malizing the MC sample to have the same data statis- 
tics at Rjj — 0, the data/MC ratio was fit to a Hnear 
function, I ± h'Rj^, where h' is the fitted slope. The re- 
gion Rjj > 1.9 was excluded because this region is more 
sensitive to energy nonlinearities and resolution. The 
Kl — > tt^tt^tt" slope parameter was then assumed to be 
hooo = h' . Applying the same procedure to our KTeV 
data yields a result consistent with the CERN-NA48 re- 
sult [5| , but with a very poor fit probability. 

In light of new information about rescattering from 
Kl — » TT+TT^TT*^, we fit our two-dimensional Dalitz plot 
to the CI3PI model [61]. With the exception of /iqoo and 
0-0 ~ 0.2, the CI3PI model parameters have been mea- 
sured or calculated theoretically, and these parameters 
are listed in Table HI In this section the fitting technique 



We first define N!^?;^^ to be the number of events recon- 
structed in a 0.05 x 0.05 Dalitz pixel (Fig. [5]) denoted by 
integers x and y. The model prediction for the number 
of events in each Dalitz pixel, N'^l'"^, is given by 



xy 



E 

y'=y-2,y + 2 



\Mooo{x', y')\'N^^,Y'S¥{x' - x,y' - y) (11) 



The quantities appearing in the prediction function are 
explained as follows. M is an overall normalization factor 
such that the total phase-space density integrals on each 
side of Eq.[TT]are the same. M.ooo{x' , y') is the matrix el- 
ement at the center of pixel x' ,y' , as calculated from the 
CI3PI model and the floated value of /looo- The remain- 
ing quantities are based on Kl — > tt^tt^tt^ MC generated 
with hooo — and no rescattering: i.e., flat phase space. 
N^y, is the number of Kl — ^ Tr'^Tr'^Tr'^ events generated in 
pixel x', y' that pass all selection criteria; note that N^y, 
is not the number of MC events reconstructed in pixel 
x' ^y' . PSF(a::' — x,y' — y) is the "pixel-spread-function," 
computed from MC, which gives the fraction of events 
generated in pixel x' ,y' that are reconstructed in pixel 
x^y. In each of the 2956 Dalitz pixels with data, the 
PSF is computed on a 5 x 5 grid around the pixel. The 
pixel size corresponds to about 5cr x 5cr in terms of the re- 
construction resolution of Xd and Yl>. On average, 70% 
of the MC events are reconstructed in the same Dalitz 
pixel as the generation pixel; 99.96% of the MC events 
are reconstructed within a 3 x 3 pixel grid centered on 
the generation pixel. 

The data are 
function 



fit with MINUIT to minimize the x 



X 



data 
xy 



xy } I xy J 



(12) 



where N^y'^'^ is the prediction function in Eq. [TI] and 
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N^^*^ is the number of reconstructed Kl 7r"7r°7r° de- 
cays in pixel x^y. The statistical uncertainty is 

Ky^r = ^.T' + (13) 

where N^p is the number of MC events reconstructed in 
pixel X, y. The two terms above represent the statistical 
uncertainty on the data and MC, respectively. 

In the fitting procedure we make an additional selec- 
tion requirement that among the three possible tt'^tt'^ 
pairings, the minimum tt^tt'^ mass, "to™o™o," is greater 
than 0.274 GeV/c^. This requirement removes 3 million 
(4.5%) Kl — > TT^TT^Tr" decays from the data sample, and 
it is applied because of a slight data-model discrepancy 
that is discussed in S lyilBl and S IVIIIAI 

The quality of the fit is illustrated by the x^- With 
/iQoo as a fit parameter and ao — 02 fixed, x^/dof = 
2805.3/2765 for all of the pixels, and Xcdgo/dof = 
125.3/130 for the subset of edge-pixels that overlap the 
Dalitz boundary. The sensitivity of Xedge illustrated 
by fitting the data without the kinematically-constrained 
energy adjustments (Eq.[Tni): in this case, x^do^c increases 
by 60. Fitting for both /looo and 00 — 02, the correspond- 
ing x^ and Xedge very similar. The results of these 
fits are presented in § IVIIIIIIXI 



VI. OBSERVATION OF INTERFERENCE 
FROM Kl TT+TT 7r° WITH RESCATTERING 

While the cusp from rescattering is clearly visible in the 
CERN-NA48 distribution of m^o^o from tt^-k^tt^ 
decays (see Fig. 2 of 0), there is no such evidence in 
our raw distribution of m™o"n from Kj^ —> tt^tt^tt^ de- 

n n n 

cays (Fig. [5)3). The rescattering effect in t: n r 

decays becomes apparent only when the data are di- 
vided by the corresponding MC distribution generated 
with pure phase-space: i.e, /iqoo — and no rescattering 
from Kl — > 7r+7r~7r° decays. These data/MC (phase- 
space) ratios are shown as data points with errors in 
Figs. \6t,d. A cusp is clearly visible in the Dalitz region 
~ 2 and m;°j"o ~ 2m^+ = 0.28 CeV. The rescatter- 
ing process tt+tt" — > tt^tt" changes from a virtual process 
(m^+TT- < 2™^^+) resulting in destructive interference, to 
a real process (m^+j^- > 2m^+) resulting in constructive 
interference. 

We use the fit results (§ IVIIip to compute a prediction 
for the data/MC (phase-space) ratio as a function of Rj^ 
and 77i"o"o ; these predictions are shown as solid curves in 
Figs.[6j:,d. The predictions agree well with our measured 
data/MC(phase-space) distributions, except for the dis- 
crepancy in the region defined by 7T7,"o"o < 0.274 GeV/c^ 
(first four bins of Fig. [5Jl). The dashed curves show 
the prediction using the CI3PI model and /iqoo replaced 
with the current PDG value, /iqoo = —0.005; these curves 
clearly do not match the KTeV distributions. To easily 
reproduce the KTeV prediction, we have parametrized 



the solid curve in Fig. [6li as a polynomial of the form: 
7^^„°^^VXo) - E ™' X - 0-281)" ' (14) 

n=0 

where ™™5Jo is the minimum tt'^tt'^ mass (GeV/c^), and 
the coefficients (m„=o,3) are given in Table [Til The root- 
mean-square precision of this parametrization is 0.023%, 
and the largest deviation of the parametrization is 0.06%. 

TABLE II: Polynomial coefficients (Eq. [14| for the 
parametrization of the solid curve in Fig.jSJi for ftooo = -1-0.59 
and m^+(ao — 02) — 0.268. Note that the coefficients depend 
on the value of m™o"o- 

7r"7r^ 



valid range of 

min tt^tt" mass mo mi m2 ma 

m™o'"o < 0.281 GeV/c^ 0.999937 3.34994 165.229 

m^o > 0.281 GeV/c^ 0.998851 0.121152 16.5534 -372.656 



VII. SYSTEMATIC UNCERTAINTIES 

Systematic uncertainties are broken into three cate- 
gories: detector & reconstruction, fitting, and external 
parameters. Within the framework of a single-parameter 
fit for hooo these categories are discussed in the subsec- 
tions below, and the systematic uncertainties on hooo are 
summarized in Table Hill The KTeV detector and analy- 
sis introduces a systematic uncertainty of 0.48 x 10~^ on 
^000- Uncertainties in external parameters, particularly 
ao — 0,2, lead to a much larger uncertainty of 1.06 x 10"'^ 
on hooo- For the two-parameter fit (/iqoo and 00 — 02; see 
§ lIXp , the systematic uncertainties are evaluated in the 
same manner, and these uncertainties are summarized 
in Table [TVl Note that when a systematic variation re- 
sults in a shift that is comparable to the the statistical 
uncertainty, we make an effort to justify an uncertainty 
that is smaller than the systematic variation; when the 
corresponding shift is much smaller than the statistical 
uncertainty, there is no need to justify a smaller uncer- 
tainty. 

A. Detector & Reconstruction 

Systematic uncertainties on hooo are mainly from ef- 
fects that bias the reconstructed Dalitz variables, Xd 
and Yd, in a manner that is not accounted for in the 
simulation. 

Kaon Scattering 

Recall that beryllium absorbers were placed 20 meters 
downstream of the primary target in order to increase 
the kaon-to-neutron ratio. Scattering in these absorbers 
affects the kaon trajectory, and hence the reconstructed 
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Fit for IIqqq (aQ-ap fixed to NA48 value) 




1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

0.5 1 1.5 2 2.5 



1 I I I I I I I I I I I I I I I I I I I I I r 

0.27 0.28 0.29 0.3 0.31 

min 7iV mass (GeV/c^) 



O 1.02 
CO 1.015 



CD 

55 0.995 
CD 

Q_ 0.99 

O 0.985 



CD 
-I—' 

CD 
Q 



0.98 



(c) 




-|- Data/MC(phase spac^) 

prediction, best-fit fiooo ^ 

prediction, h^gg = -0.005 



I I I I I I I I I I I I I I I I I I I I I I I I 
0.5 1 1.5 2 2.5 



1.02 
1.015 



0.98 



(d) 




fv -\- Data/IVIC(plnase space) 

y prediction, best-fit h^^^ 

prediction, h^^^ = -0.005 



I I I I I I I I I I I — I I I I I I I I I I 
0.27 0.28 0.29 0.3 0.31 

min TiV mass (GeV/c^) 



FIG. 6; For the 68.3 million Kl — > Tr^Tr^Tr" in the KTeV sample, projected Dalitz distributions are shown for (a) and (b) 
m™o"o- The average reconstruction resolution determined by the simulation is a{R^) ~ 0.014 and a{mm m^o^o) ~ 0.3 MeV/c^: 
these resolutions are indicated by a IOct marker on each plot. The data/MC(phase-space) ratio is shown as a function of (c) 
Ro and (d) m™o"o (points with error bars). The solid curve is the prediction from our best fit hooo- The dashed curve is the 
prediction using ftooo (PDG06) = (—5. Oil. 4) x 10~^. The arrow in (d) shows the selection requirement m™o"o > 0.274 GeV/c^. 
Note that previous analyses [j, 01 ignored rescattering and excluded R^ > 1.9 ; the corresponding data/MC ratio was assumed 
to be a straight line with slope of —0.005. 



Dalitz variables. If absorber-scattering is turned off in 
the MC, the resulting value of /iqoo changes by 0.5 x 
10"'^. Based on studies of kaon trajectories with K'^ 
7r"*'7r~ decays in the vacuum beam, we assign a systematic 



uncertainty on /iqoo equal to 10% of the change when 
scattering is turned off in the simulation: 0.05 x 10~^. 

Accidental Activity 

Energy deposits from accidental activity in the Csl 
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calorimeter can modify the reconstructed photon ener- 
gies. In the reconstruction, events are rejected if any 
of the six photon clusters has accidental activity within 
a 19 nanosecond window prior to the start-time of the 
event. Removing this cut increases the level of acciden- 
tal activity, and changes /looo by 0.02 x 10^"^; we include 
this difference as a systematic uncertainty. 

Photon Energy Scale 

The photon energy scale is determined in the i?e(e'/e) 
analysis by comparing the data and MC vertex distribu- 
tions for K'^ tt'^tt^ decays downstream of the regener- 
ator. These decays are mainly due to the i^s-component 
of the neutral kaon. The active veto system rejects de- 
cays inside the regenerator, resulting in a rapidly rising 
distribution just downstream of the regenerator. The 
data-MC vertex comparison has a discrepancy that is 
slightly dependent on kaon energy, and the magnitude of 
the discrepancy no more than 3 cm; this data-MC shift in 
the vertex corresponds to an energy-scale discrepancy of 
up to ^ 0.05%. An energy scale correction is empirically 
derived to remove this small discrepancy in — > tt'^tt'^ 
decays, and this "tt^tt^" correction is applied to photon 
energies in the Kl tt^tt^tt'^ Dalitz analysis. As a sys- 
tematic test, the Dalitz analysis is performed with no 
energy scale correction: /iqoo changes by 0.06 x 10^'^ and 
is included as a systematic error. 



TABLE III: Systematic uncertainties on /iqoo- For each ex- 
ternal parameter X, the sign (-1- or — ) is indicated for the 
partial derivative, dhooo/dX , so that our /iqoo result can be 
updated when an external parameter is updated. 



source of 


uncertainty on /looo 


uncertainty 


(xlQ-^) 


DETECTOR & RECON 




kaon scattering 


0.05 


accidentals 


0.02 


photon energy scale 


0.06 


energy resolution 


0.04 


low-side energy tail 


0.02 


position resolution 


0.07 


xl-cut 


0.07 


(sub-total) 


(0.13) 


FITTING 




MC statistics 


0.14 


Ignore PSF for NEy""^ 


0.02 


remove m™d"o cut 


0.44 


(sub-total) 


(0.46) 


KTeV TOTAL 


0.48 




EXTERNAL 




(ao - a2)m^+ 


(+) 1.03 




(-) 0.12 


ro, r2 


(+) 0.21, (+) 0.04 


At/Al 


(+) 0.01 


g+-o, h+-o 


(-) 0.05, (-) 0.05 


(sub-total) 


(1.06) 



Photon Energy Resolution 

The simulated energy resolution is adjusted by about 
0.3% to match the energy resolution for electrons from 
Kl TT^e^v decays. The resulting photon energy 
resolution is well simulated, as illustrated by the excel- 
lent data-MC agreement in the 7r°7r°7r'^-mass distribution 
(Fig. UK). As a systematic test, we increase the simulated 
resolution by an additional 0.3%: the change in /iqoo is 
0.04 X 10"^, and is included as a systematic uncertainty. 

Low-Side Energy Tail 

The effects of photo-nuclear interactions and wrapping 
material in the Csl calorimeter can result in photon en- 
ergies measured well below a few-sigma fluctuation in the 
expected photostatistics. As described in § IIIII this non- 
Gaussian tail has been measured using electrons from 
Kl TT^e^v decays, and modeled in the simulation. 
Based on the data-MC agreement in the low-side E/p tail 
for electrons, we assign a 20% uncertainty on our under- 
standing of this effect. As an illustration, note that the 
Gaussian energy resolution (0.8%) predicts that 0.02% 
of the photons will be reconstructed with an energy that 
is at least 3% below the true value; the effect of the 
non-Gaussian tail is that 0.8% of the reconstructed pho- 
ton energies are at least 3% low. As a systematic test, 
we remove simulated decays in which any photon loses 
more than 3% of its energy due to this non-Gaussian 
process. This test rejects 6 x 0.8% ~ 5% of the generated 
Kl TT^TT^TT^ decays. After applying selection require- 
ments, the MC sample is reduced by 2%, which is smaller 
than the reduction for generated decays. Using this test- 
MC sample, the change in /looo is 0.1 x 10~^ compared 
to using the nominal MC; as explained above, we include 
20% of this change, 0.02 x lO^'^, as a systematic uncer- 
tainty on /looo- 

Photon Position Resolution 

Turning off the "un-smearing" (§ IIII[) of the MC photon 
positions results in a change of 0.3 x 10^'' in /iqoo- Based 
on the data-MC agreement in the electron position res- 
olution from Kl n^e^v decays, we take 20% of this 
change, 0.07 x 10^'^, as a systematic uncertainty. 

Xl Cut ^ ^ 

The determination of the Dalitz variables is performed 
using adjusted photon energies, where the adjustment is 
done for each Kl —> tt'^tt^tt^ decay by minimizing the 
"energy-xl;" in Eq. [TOl The selection requirement is 
Xe < 50. As a systematic test, this cut is relaxed to 
X% < 1000; the change in /iqoo is 0.07 x 10~^, and is 
included as a systematic uncertainty. 

B. Fitting 



MC Statistics 

The simulated sample consists of 124.9 million Kl 
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7r°7r°7r° decays that satisfy the selection requirements 
(1.8x the data statistics). This sample results in a MC- 
statistics uncertainty of 0.14 x lO^'^ on ft-ooo- 

Pixel Migration 

The reconstructed pixel location in the Dalitz plot 
(XotYd) can be different than the true pixel location 
This pixel migration is accounted for by using the pixel- 
spread-function (PSF) in Eq. [TT] to predict the number 
of reconstructed K]^ tt^tt^tt^ decays in each pixel. As 
a systematic test, we ignore pixel migration by setting 
PSF(a;' ~x,y'-y) = 6{x' -x,y' - y) and replacing N^^, 
with the number of events reconstructed in each pixel; 
the change in /loooi 0.02 x 10"'^, is included as a system- 
atic uncertainty. 

Data-Model Discrepancy 

As shown in Fig.[6]i, the Dalitz region defined by "^^"o < 
0.274 GeV/c^ shows a data-model discrepancy, and this 
region is therefore excluded from the nominal fit. Includ- 
ing this region in the fit changes ft-ooo by 0.44 x 10^"^, 
and we include this difference as a systematic uncer- 
tainty. Additional discussion on this discrepancy is given 
in S IVIII Al 



C. External Parameters 

The CI3PI model depends on several parameters listed 
in Table [H The uncertainties in these parameters have 
been propagated through the /iqoo fit- The net /iqoo un- 
certainty from these external parameters is 1.06 x 10^"^. 
This uncertainty is almost entirely due to the uncertainty 
in the difference in scattering lengths, 09 — 02. 

In Table IIIIl we have also included the sign of each 
partial derivative so that our /looo result can be updated 
when an external parameter is updated. For example, 
c?^ooo/c?(ao"^7r+) = —0.12/0.013, where the numerator 
and denominator are from Tables Hill and HI respectively. 



VIII. RESULT FOR ftooo WITH FIXED ao - aa 

Here we fix 111^+ (00—02) = 0.268 as measured by NA48 
Q), and determine /iqoo- The result from minimizing the 
in Eq.IIlis 

/looo = (+0.59 ± 0.20,ta0 X 10"^ (15) 
xVdof = 2805.3/2765 (all pixels) (16) 
xVdof = 125.3/130 (edge pixels) . (17) 

where the statistical uncertainty is from 68.3 million de- 
cays in the data sample. To check our modeling near the 
Dalitz boundary, the x^ is shown for the subset of "edge 
pixels" that overlap the Dalitz boundary. 



Including the systematic uncertainty, the final result 
for the quadratic slope parameter is 

hmo = (+0.59 ± 0.20,tat ± 0.48,^,,* ± 1.06^.*) x lO^^ 

(18) 

= (+0.59+1.19) X 10"^ (19) 

where the uncertainties are from data statistics, KTeV 
systematic errors, and external systematics errors. 

A. Crosschecks on /looo 

Some crosschecks on the result for /looo a-re shown 
in Fig. [71 The separate measurements for each year 
are consistent, as well as the separate measurements 
from each vacuum beam. The last crosscheck involves 
the asymmetry between the minimum and maximum 
photon energy, which could expose potential problems 
related to non-linearities in the photon energy measure- 
ment. The ratio between the minimum and maximum 
photon energies, = E"^^'^ / E'^^^ , is used to define 
five sub-samples with roughly equal statistics: r-y = 
{0, 0.14}, {0.14, 0.18}, {0.18, 0.23}, {0.23, 0.31}, {0.31, 1}. 
The five independent measurements of /iqoo are consis- 
tent. 

Concerning the data-model discrepancy in the Dalitz 
plot region m^j^ < 0.274 GeV/c^ (Fig. Hi), we have 
performed many checks to investigate if the problem is 
related to our analysis. For example, the MC energy res- 
olution was degraded by an additional 0.8%, an extreme 
change that is nearly three times larger than the stan- 
dard 0.3% smearing: the corresponding change in /iqoo is 
1.6 times the statistical uncertainty {(y^tat)^ but the data- 
model discrepancy remains unchanged. In another test, 
an extreme energy nonlinearity of 0.3% per 100 GeV is 
introduced into the simulated energy measurements; 
changes by 0.5cr^(^j, and the data-model discrepancy is 
again unchanged. These highly exaggerated tests suggest 
that the KTeV energy reconstruction is not responsible 
for the data-model discrepancy. We have also checked 
that the data-model discrepancy is unchanged for the 
following tests: vary best Xpair cut between 4 and 100 
(nominal cut is 10), remove requirement that the sec- 
ond smallest Xpair value is greater than 30, allow no hits 
and up to six hits in the scintillator hodoscope (to check 
photon conversions), allow photons to hit a Csl crystal 
adjacent to the beam holes (Fig. |3l), remove requirement 
on Csl cluster energy deposited before event (increases 
effect from accidentals), vary cut on Xe frorn < 10 to 
no cut (Fig. Sb), remove simulated decays in which any 
photon loses more than 3% of its energy in the Csl (see 
systematic test "Low-Side Energy Tail" in § IVII Ap , use 
reconstructed Csl photon energies instead of adjusted en- 
ergies based on kinematic constraints. 

Photon conversions in detector material result in e^e^ 
pairs that are reconstructed as a single photon. A scin- 
tillator hodoscope just upstream of the Csl calorimeter 
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FIG. 7: Crosscheck measurements of hooo'- data-taking years 
(96,97,99), left and right vacuum beams (L,R), and min/max 
photon-energy ratio as discussed in the text. 

Measurements within each category (between vertical lines) 
are statistically independent. Error bars reflect the statistical 
uncertainties from the data and MC samples. 



tags such e+e^ pah's. The standard analysis allows up 
to one hit in this hodoscope. As a systematic test, we 
compare results with (i) no requirement on hodoscope 
hits, and with (ii) a requirement that there are no hits 
in the hodoscope. For these two samples, there is a 15% 
difference in the number of reconstructed Kj^ — > tt^tt'^tt'^ 
decays, and the difference in /iqoo is (0.07 ±0.09) x 10~^. 

As a final crosscheck, the analysis is repeated using the 
reconstructed photon energies instead of the adjusted en- 
ergies based on kinematic constraints from the and 
7r° masses (see Xe Eq.fTOl). Using unconstrained Dalitz 
variables, the resulting value of /iqoo changes by 1.2a'^^^^ 
compared to the nominal result. However, compared to 
the nominal result in Eqs. [TBlfTTl the overall fit-x^ in- 
creases by 120, and the fit-x^ for the edge pixels increases 
by nearly 60. This increase in indicates that the reso- 
lution is not modeled as well for the unconstrained Dalitz 
variables, and it illustrates the importance of the kine- 
matic constraints. 



IX. MEASUREMENT OF ao - aa AND /looo 
WITH Kl TT W DECAYS 

Here we use Kl — > tt^tt^tt" decays to measure both 
the quadratic slope parameter and the difference in pion 
scattering lengths. The fit procedure is described in §|Vl 
but now we float oq — 02 instead of fixing it to the value 
measured by NA48 0. Fitting our data for both ft-ooo 
and ao — 02 in a two-parameter fit, we find 

haoa = 

{-2m±0.62stat ± 0.72,j,,t ± 0.28e.t) X 10-3 (20) 
= (-2.09 ±0.99) X 10-3 (21) 



m^+{ao-a2) = 
0.215 ±0.014,tat ± 0.025,j,,t ± 0.006ext (22) 
= 0.215 ±0.031 (23) 

Pha = +0.939 (24) 
xVdof = 2790.6/2764 (all pixels) (25) 
xVdof = 126.3/130 (edge pixels) . (26) 

The uncertainties are from data statistics, KTeV sys- 
tematic errors, and external systematic errors. The sys- 
tematic uncertainties are evaluated in the same manner 
as for the one-parameter fit for /iqoo (§ IVH[) : these un- 
certainties are summarized in Table ITVl The data-model 
comparisons are shown in Fig. [51 

Compared to the fit in which ao — 02 is fixed (Eq. llSp . 
the statistical uncertainty on /looo is more than x 3 larger 
but the total uncertainty is slightly smaller. The reason 
for the smaller /looo-uncertainty when ao — 02 is floated 
is related to the nonlinear dependence of the correlation 
between /looo and ao — a2. When ao — 02 = 0.268 is 
fixed , dhoon/d{af) — 02) ~ 0.06. For our best-fit value of 
flo — a2 = 0.215, dhoQo/d{ao — 02) — 0.04 and hence ft-000 
is less sensitive to variations in gq — a2. The asymmetry 
between +lcr and —la variations is about 10%, so we 
simply averaged the ±lcr variations and quote symmetric 
errors. 



A. Comparisons of Results 

We begin by comparing the hooo result for the two 
different fits. Compared to the one-parameter fit where 
ao — a2 is fixed, the statistical uncertainty on hooo from 
the two-parameter fit (Eq. [^0)1 is about x3 larger and 
the KTeV systematic uncertainty is xl.5 larger. The 
systematic uncertainty increases by less than the sta- 
tistical uncertainty because the largest source of uncer- 
tainty (cut on m™o'"o) is similar in both the one- and 
two-parameter fits. While the ft-ooo measurement errors 
are much larger for the two-parameter fit, the external 
uncertainty is x4 smaller than the external uncertainty 
for the one-parameter fit. The large difference in the 
external uncertainties is driven by the large correlation 
(pha — +0.939) between hooo and ao — a2. The over- 
all uncertainty on hooo is nearly the same for the one- 
and two-parameter fits; after accounting for the different 
sources of uncertainty in each fit, the significance on the 
different values of hooo (+0.59 vs. —2.09) is estimated to 
be 2(7. 

Next we compare our ao — a2 result to the NA48 anal- 
ysis based on Tr^Tr'^Tr" decays where they reported 
ag - 02 = 0.268 ± 0.017. The KTeV statistical uncer- 
tainty on ao — a2 is about 40% larger [l^ even though our 
Kl tt^tt^tt" sample is more than twice as large as their 
(NA48) K^ 7r*7r°7r° sample; the larger statistical un- 
certainty from Kl tt^tt^tt^ decays is due to the much 
smaller rescattering effect compared to K^ — > tt^tt^tt^ 
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decays. Our overall uncertainty on gq — ai is nearly x2 
larger than that obtained by NA48. The KTeV and 
NA48 results on — are consistent at the level of 
1.5(7. Our result is also compatible with the DIRAC re- 
sult based on measurin g th e lifetime of the tt+tt^ atom: 
|ao~a2| = 0.264lE]°33g|. 

TABLE IV: Systematic uncertainties on /iqoo and (ao — 
a2). For each external parameter X, the sign (+ or — ) is 
indicated for the partial derivative, dh^Q^I dX , so that our 



results can be updated if an 


external paramet 


er is updated. 


source of 


uncertainty on 


uncertainty 


10'' X /looo 


m„+{ao - 02) 


DETECTOR & RECON 






kaon scattering 


0.05 


0.000 


accidentals 


0.03 


0.000 


photon energy scale 


0.13 


0.002 


energy resolution 


0.11 


0.002 


low-side energy tail 


0.02 


0.000 


position resolution 


0.10 


0.001 


xl-cut 


0.07 


0.000 


(sub-total) 


(0.21) 


(0.003) 


FITTING 






MC statistics 


0.46 


0.011 


Ignore PSF for N'^-^'^ 


0.07 


0.002 


remove m™o"o cut 


0.51 


0.022 


KTeV TOTAL 


0.72 


0.025 




EXTERNAL 








(-) 0.07 


(-) 0.002 


tq 


(+) 0.15 


(+) 0.001 


ri 


(-) 0.01 


(-) 0.000 


At/Al 


(-) 0.08 


(-) 0.002 


g+-o 


{+) 0.12 


{+) 0.004 




(-) 0.17 


(-) 0.003 


(sub-total) 


(0.28) 


(0.006) 



X. CONCLUSION 

We have made the first observation of interference be- 
tween the Kl —>■ TT^Tr'^Tr'' decay amplitude, and the ampli- 
tude for Kl 7r+7r~7r° with the final-state rescattering 



process tt+tt" tt^tt". When comparing our data to a 
Monte Carlo sample of Kl Tr^Tr^Tr" decays generated 
with pure phase-space, we see a cusp in the data/MC 
distribution-ratio of minimum tt'^tt'^ mass. This cusp is 
not visible in the data distribution (Fig. [61d); rather, it is 
visible only in the data/MC ratio (Fig.[6ji). 

Using the CI3PI model [6^ to account for rescatter- 
ing, and fixing 00 — 02 to the value measured with — > 
7r^7r°7r° decays 0, we have measured the quadratic slope 
parameter, /iqoo = (-1-0.59 ± 1.19) x 10"^, where the 
largest source of uncertainty is from the uncertainty on 
09 — 02. This result is consistent with zero, and it dis- 
agrees with the average of previous measurements that 
did not account for rescattering. The CI3PI model de- 
scribes the data well for most of the Kl Tr^Tr^Tr" 
phase space, but there is a notable 0.3% discrepancy in 
the region where the minimum 7r°7r° mass is less than 
0.274 GeV/c^. We have excluded this discrepant region 
from our nominal fits, but have included this region to 
evaluate systematic uncertainties. To investigate the pos- 
sibility that the data-model discrepancy is from our anal- 
ysis, we have made extreme variations in the simulation 
of the photon energy scale and resolution (§ IVIII Ap and 
found that such drastic changes have no impact on the 
discrepancy. We have not been able to numerically verify 
the calculation of the model, but for future comparisons 
we have left a convenient parametrization (Eq. [TJ] and 
Table ini). 

We have repeated our phase space analysis by floating 
oq — 02 rather than fixing it to the value reported by 
NA48. Detailed results are presented in § |lXl Our value 
of Og — 02 is consistent with that found by NA48, but 
with an uncertainty that is nearly twice as large. 
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Fit for h^^^ and a^-a 
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FIG. 8: Same as Fig. |6j:-d, except ao — 02 is floated in the fit instead of fixed to the NA48 value. 



